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PREFACE 

This  report  is  an  integration  ;  d  extension  of  concepts  developed  under  several  earlier  Visibility 
Laboratory  contracts. 

The  first  of  these,  N5ori-07831,  N5ori-07864,  and  NObs-50378  were  contracts  between  the  United 
States  Navy  and  the  Massachusetts  Institute  of  Technology  when,  prior  to  1952,  the  Visibility  Labor¬ 
atory  was  a  part  of  that  institution  (see  reference  1,  entitled  "The  Visibility  of  Submerged  Objects”). 
Another,  much  later  contract,  NObs-34075  (Assignment  1),  concerned  the  night  aerjal  search  for  avi¬ 
ators  downed  at  sea  for  the  U.  S.  Navy  Bureau  of  Naval  Weapons.  Three  other  contracts,  NObs-84075 
(Assignment  3),  NObs-92058  (Task  1),  and  NObs-95251  (Task  1)  concerned  ocean  surveillance  from 
satellites  for  the  Bureau  of  Ships.  Still  another,  Nonr-2216  (Tasks  25  and  27),  studied  the  appearance 
of  water  wave  patterns  for  the  Office  of  Naval  Research.  Finally  contract  NObs-84075  (Assignment 
9)  provided  for  the  development  and  reporting  of  technical  materials  useful  for  visibility  calculations. 

The  writing  of  this  report  was  accomplished  partially  under  Assignment  9  of  Contract  NObs-84075 
and  partially  under  Task  1  of  NObs-95251 .  Final  editing,  reproduction  and  distribution  was  accom¬ 
plished  under  Contract  N00024-68-C-1100,  Task  1  to  which  the  subject  matter  of  the  report  contributes 
directly  as  an  integral  part  of  the  Laboratory's  "investigation  of  the  detectability  and  recognizability 
of  military  targets  in  the  sea  environment". 
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DIRECTIONAL  RADIANCE 
(LUMINANCE)  OF  THE  SEA  SURFACE 

Jacqueline  I.  Gordon 


1.  INTRODUCTION 

Equations  for  the  computation  of  directional  sea  radiance  and/or  luminance  are  presented.  Three  types 
of  conditions  for  which  the  calculations  may  be  made  are  each  treated  separately:  first,  the  radiance  or 
luminance  of  the  calm  sea:  second,  the  unresolved  radiance  or  luminance  of  the  time-  or  space-averaged 
wind-ruffled  sea;  and  third,  the  radiance  or  luminance  of  the  sea  surface  in  which  the  large  wave  pattern  _ 
is  resolved  but  the  small  capillary  waves  are  unresolved. 

This  report  is  for  the  purpose  of  presenting  the  complete  set  of  equations,  with  derivations,  that  are 
needed  for  the  calculation  of  the  directional  sea  radiance  (luminance).  The  basis  for  the  calculations  in 
Sections  2  and  3  can  be  found  in  an  earlier  report  by  Duntley1.  The  extensions  and  modifications  are  in 
agreement  with  the  theories  and  equations  given  by  Preisendorfer2,  but  the  treatment  of  the  details  of  the 
sea  radiance  (luminance)  computations  given  herein  are  in  greater  detail.  Many  of  the  equations  in  Section 
4  also  appear  in  the  unclassified  section  of  a  final  report  for  an  ONR  project3;  no  derivations  were  given 
there,  however,  so  the  present  treatment,  while  parallel,  contains  modifications  and  is  more  complete. 

Angular  notation  throughout  is  as  follows:  The  zenith  angle  is  designated  by  6.  See  Fig.  1.  The  azi¬ 
muth  is  with  reference  to  the  azimuth  of  the  sun  and  is  noted  by  <A  These  angles  are  always  given  as 
parenthetic  attachments  to  a  radiance  N  or  a  luminance  B,  the  angles  thus  denoting  the  direction  of  the 
path  of  sight  of  the  photometer.  The  altitude  is  not  included  in  the  parenthesis  since  it  is  always  sea 
level  or  zero  in  an  inherent  sea  radiance  or  luminance  computation 


Radiance  units  used  throughout  are  watts/Q  sq.  ft.;  irradiance  units  are  watts/sq.  ft.  Luminance  units 
are  in  Ft-Lamberts  and  illuminance  units  are  in  lumens/sq.  ft.  The  directional  radiance  of  the  sea  surface 
would  thus  be  N(0p,<£p)  and  the  luminance  B{6p,<f>p). 


Zenith 


Fig.  1.  Plane  passed  through  path  of  sight  and  zenith. 


The  inherent  radiance  (luminance)  of  calm  water  will  be  discussed  first,  followed  by  the  radiance 
(luminance)  of  the  time-  or  space-averaged  wind-ruffled  sea  surface.  Next  will  be  described  the  exten¬ 
sion  of  these  ideas  to  approximate  the  radiance  (luminance)  of  the  sea  surface  as  it  might  appear  in  a 
photograph  with,  for  example,  a  resolution  of  1  ft. 

A  program  for  the  CDC  3600  computer  has  been  completed  which  carries  out  the  computations  herein 
described  for  calm  water,  the  time-averaged  wind-ruffled  sea,  and  the  tilted  sea  surface  of  intermediate  re¬ 
solution  (in  time-  or  space-averaging).  This  report  contains  the  equations  which  are  the  basis  of 
that  program. 


2.  CALM  WATER 

The  directional  sea  radiance  (luminance)  for  calm  water  is  the  sum  of  the  reflected  radiance 
(luminance)  of  the  sun  or  the  sky  and  the  upwelling  radiance  (luminance).  The  radiance  upwelling  from 
within  the  sea  results  from  sunlight  and  skylight  which  has  been  transmitted  downward  through  the  air- 
water  surface,  scattered  upward  by  the  (optically  infinite)  deep  water,  and  then  I  ansmitted  through  the 
water-air  surface  toward  the  sensor. 

The  reflectance  and  transmittance  of  the  air-water  surface  will  be  reviewed  first,  followed  by  a  dis¬ 
cussion  of  the  index  of  refraction  of  the  ocean  surface.  Next,  the  equations  for  computing  the  sky  or  sun 
reflectance  will  be  presented.  Then  the  upwelling  radiance  (luminance)  will  be  investigated. 

Most  of  the  equations  in  this  section  come  directly  from  Reference  1.  They  are  reproduced  here  not 
only  for  convenience  but  also  to  allow  certain  necessary  modifications  to  be  made. 

2.1  Fresnel  Reflectance  and  Transmittance 

The  Fresnel  reflectance,  r(0),  and  transmittance,  t (0),  are  dependent  on  the  index  of  refraction  of  the 
water,  n1,  the  index  of  refraction  of  the  air,  n,  and  the  zenith  angle  of  the  incident  radiation,  0. 

The  refracted  angle,  0,  see  Fig.  2,  beneath  the  water  surface  is  defined  by  Snell's  Law  to  be 


2 


0  =  sin  1  [(n/n1 }  sin  0] ' 


(2.1.1) 


Fig.  2.  Reflection,  transmission,  and  refraction. 


The  Fresnel  reflectance,  r (0),  can  be  found  as  follows: 


r(0)  =  [tan2  (0  -  0)/2  tan2  (0  +  0))  +  [sin2  (0  -  0)/2  sin2  (0  +  0)]  (2.1.2) 


The  proportion  of  the  radiation  transmitted  through  the  surface  is. 


t(0)  =  1  -  r(0) 


(2.1.3) 


The  index  of  refraction  of  air,  n  =  1. 


2.2  Index  of  Refraction  of  Sea  Water 

The  index  of  refraction  of  the  surface  sea  water  depends  upon  the  wavelength  of  the  incident  radiation; 
and  upon  the  water  temperature  and  salinity.  The  temperature  range  is  -2C  to  +30°C  The  normal  salinity 
range  is  from  33  to  37  parts  per  mi  I  le. 4  The  effect  of  salinity  and  temperature  in  these  ranges  is  an  order 
of  magnitude  less  than  the  effect  of  wavelength  of  the  incident  light.  Therefore  it  is  reasonable  to  use  the 
index  of  refraction  for  water  with  a  salinity  of  35  parts  per  mi  He  which  is  the  average  for  the  open  ocean 
Also  it  is  reasonable  to  use  a  water  temperature  of  20rC  (only  other  data  by  wavelength  available  was  at 
0r  and  40'  C) 


The  index  of  refraction  at  various  wavelengths  is  given  in  Table  1  for  distilled  water  at  20°C.S 
The  index  of  refn~tion  for  sea  water  with  3  iity  of  35  parts  per  mi  tie  at  20°C  is  given  in  the  same 
table  for  the  sodium-0  wavelength,  0.589  ftm.6  The  values  for  the  other  wavelengths  for  sea  water  were 
estimated  by  multiplying  the  distilled  water  values  by  factors.  A  factor  was  obtained  for  each  wave¬ 
length  from  the  distilled  water  data  at  20°C  from  the  ratio  of  the  index  for  that  wavelength  to  the  index 
for  sodium.  The  factor  for  0.546  pm,  mercury,  was  given  in  Ref.  6 

The  indices  in  Table  1  are  depicted  graphically  in  Fig.  3.  Indices  of  refraction  for  intermediate  wave¬ 
lengths  can  be  found  by  interpolating  between  the  values  given.  The  suggested  index  of  refraction  for 
photopic  calculations  for  seawater  is  1.341. 


Table  1.  Indices  of  refraction  of  water.* 
Temperature  -Salinity  Wavelength  m) 


°C 

°/oo 

0.397 

0.434 

0.486 

0.546 

0.589 

0.656 

Distilled  water 

20 

0 

1.3435 

1.3404 

1.3372 

1.3345 

1.3330 

1.3312 

Sea  water 

20 

35 

(1.3500) 

(1.3468) 

(1.3436) 

(1.3409) 

1.3394 

(1.3376) 

Values  in  parentheses  are  estimated. 


Fig.  3.  Index  of  refraction  of  sea  water  and  distilled  water. 


2.3  Sky  and  Sun  Reflectance 

The  calm  water  reflects  the  radiance  from  the  sky  or  sun  which  is  at  the  angular  position  0B,  <£Binto 
the  path  of  sight  at  0  ,  $p. 


sky  o 


r  tun  (0p.  <£P>  =  N(0b.  <£b>  r(0B) 


(2.3.1) 


The  reflectance  r(0B)  is  the  Fresnel  reflectance  at  angle  0B. 

The  sky  or  sun  position  0D,  <f>B  is  in  the  same  azimuth  as  path  of  sight  $p  but  at  the  zenith  angle  of. 
the  observer  position  0o.  The  path  of  sight  zenith  angle,  6p,  and  azimuth,  $p,  are  in  the  opposite  direc¬ 
tion  from  the  sensor  or  observer  position  zenith  angle,  6o,  and  azimuth,  <t>0 ■  Thus, 


» 


<f>B  =  <f>P  =^o  +  180°  (2.3.2) 

0B  =  0O  =  180° -0p  (2.3.3) 


Fig.  4.  Plane  passed  through  the  path  of  sight  and  the  zenith. 


Equation  (2.3.1)  can  be  rewritten  in  terms  of  0o  and  <f>p. 


N 


sky 


o,  .u,  UV  <V  =  N {0o.  6'  r(0j 


(2.3.4) 
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(2.3.5) 


The  sun  will  be  reflected)  when  6o  =  0s  and  0p  -  0;  thus 

N...  O8OO-0..  0)  =  Ns,  (0..O)r(0,) 
where  Nt  (0.,  0)  is  the  apparent  sun  radiance  averaged  over  the  solar  disk. 

The  photopic  equations  are  directly  analogous  as  follows: 

B.xy  .,™„(*p.*p)  =  B(0o.*p)  r(0o)  (2.3.6) 

and 

B.un  (180°  -  e%,  0)  =  B,  (0.,  0)  r(0.)  (2.3.7) 


2.4  Upwelling  Radiance  (Luminance) 

The  upwelling  radiance  (luminance)  is  light  from  the  sun  and  sky  which  has  gone  through  three  pro¬ 
cesses.  First,  the  light  is  transmitted  through  the  air-water  surface.  The  downwelling  irradiance  (illumi¬ 
nance)  is  the  light  measured  just  beneath  the  ocean  surface  by  an  irradiometer  (illuminometer).  Second, 
some  of  the  downwelling  irradiance  is  reflected  by  the  (optically  infinite)  deepwater,  and  interreflected 
between  the  water-air  surface  and  the  deep  water.  The  upwelling  radiance  (luminance)  just  beneath  the 
water  surface  in  the  appropriate  direction  for  the  sensor's  (observer's)  path  of  sight,  results  from  these  re¬ 
flections.  Third,  the  upwelling  radiance  (luminance)  is  transmitted  through  the  water-air  surface  into  the 
path  of  sight  of  the  sensor  (observer).  This  is  the  upwelling  radiance  (luminance)  as  measured  just  above 
the  water  surface. 


2.4.1  DOWNWELLING  IRRADIANCE  (ILLUMINANCE)  BENEATH  THE  WATER  SURFACE. 


The  sky  radiance,  N(0b,i/>b),  from  a  solid  angle,  dO,  and  projected  area,  ds  cos0B,  enters  the  water  at 
an  angle  0B.  The  proportion  entering  is  t(0B). 


Fig.  5.  Transmittance  geometry. 
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The  radiance  just  beneath  the  surface,  N(0B,  <£B)  is  related  to  the  radiance  above  the  surface  by  the  ratio 
of  the  respective  solid  angle*  and  projected  areas,  and  by  the  transmittance  t(0B),  see  Fig.  5. 


N(0b,  <f>B)  =  N(0B,  0b)  t(0B) 


ds  cos0B  d!2 
ds  cos0B  dQ’ 


(2.4.1) 


A  flat  plate  collector  will  receive  this  sky  radiance,  N(0B,  <£B),  at  the  refracted  angle,  0B.  The  col¬ 
lector  will  weight  the  amount  received  by  cos0BdO'.  The  downwelling  irradiance  beneath  the  water  sur¬ 


face,  HD(0, 0)  is  thus, 


HD(0,0)  =  /„ N(®b-  <£B)cos0BdQ' 


(2.4.2) 


Substituting  from  Eq.  (2.4.1)  for  N(@B, 


HD(0,0)  /2^N(0b,  0B)t(0B)cos0BdO 


(2.4.3) 


The  downwelling  irradiance  can  also  be  expressed  as  a  double  summation,  with  AflB  =  sin0BA0BA<£B. 
Thus,  Eq.  2.4.3  becomes. 


Hd(0,0)  =  N(0B1,  <£Bj)  A<£BJ]  t(0B1)  cos0m  s\nem  MBl  (2.4.4) 


The  sun  contribution  can  be  treated  as  a  separate  term,  in  an  analogous  fashion.  Therefore,  the  total 
downwelling  irradiance  from  sun  and  sky  would  be, 

hd(°-  0)  =  jw N(0b,  0b)  t(0B)  cos0B  dfiB  +  NS(0S,  0)  t(0s)  cos0s  das  (2.4.5) 

where  dftt  is  the  solid  angle  of  the  solar  disk;  dfi,  =  6.8193  x  10-5  steradians  at  mean  solar  distance. 
The  maximum  variability  in  dfit,  at  the  perihelion  and  aphelion,  is  +  3.4  percent. 

The  equation  for  the  downwelling  illuminance  beneath  the  water  surface  has  an  analogous  derivation; 
but  it  contains  a  n  term  since  the  luminances  used  are  in  ft-L. 


Ed(0.0) 


■% 


2n  B(0b,  0b>  cos0BdfiB  +  BS(0S,O)  t(0s)  cos0s  d«s 


(2.4.6) 


2.4.2  UPWELUNG  RADIANCE  (LUMINANCE)  BENEATH  THE  WATER  SURFACE 

The  reflectance  of  the  (optically  infinite)  deep  water  is  a  function  of  wavelength.  It  is  measured  just 
beneath  the  water  surface  by  two  flat  plate  collectors:  one  measuring  the  upwelling  spectral  irradiance 
^(180,0)  ana  one  measuring  the  downwelling  spectral  irradiance  HD(0,0). 


*  The  parenthetic  expression  (0,0)  denotes  the  zenith  angle  and  azimuth  of  the  normal  vector  from  the  surface  of  the 
irrndionietor. 
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(2.4.7) 


^  (180,0) 
Hd(0.0) 


The  photopic  deep  water  reflectance  is  defined  by  the  ratio  of  the  upwelling  illuminance  to  the  down- 


welling  illuminance. 


E„  (180.0) 
ED(0.0) 


(2.4.8) 


The  deep  water  reflectance  is  a  function  of  the  water  clarity  as  well  as  wavelength.  The  standard  photopic 
value  for  deep  water  reflectance  is  0.02  but  values  ranging  from  0.01  to  0.15  have  been  reported.7 

The  reflectance  of  deep  water  as  a  function  of  wavelength  is  given  in  Table  2,  and  Fig.  6.  The  curve 
represents  measurements  made  with  a  spectrogeograph  in  a  glass-bottomed  boat  in  the  Gulf  Stream  off 
Dania,  Florida. *-9-10  Values  below  0.420/<m  and  above  0.570 /xm  are  extrapolated.  Note  that 
RJ0.550)  ^  0.02.  The  photopic  deep  water  reflectance  of  0.02  can  be  computed  from  the  following  equa¬ 
tion;  where  the  incoming  radiation  is  the  solar  spectral  radiation  from  Johnson. 11 

4  R~WHAVAdA 

Roo(y>  =  —7s -  (2-4.9) 

J0  HA  Va  dA 

The  reflectance  of  (optically  infinite)  deep  water  for  a  radiometric  sensor  having  a  broad  band  spectral 
sensitivity  of  would  be  computed  in  a  similar  manner. 


Table  2.  Reflectance  of  (optically  infinite)  deep  water.* 


Wavelength 

(//m) 

■u*> 

Wavelength 

(ftm) 

Wavelength 

(/xm) 

0.40 

(0.075) 

0.50 

0.053 

0.60 

(0.0035) 

0.41 

(0.079) 

0.51 

0.040 

0.61 

(0.0029) 

0.42 

0.082 

0.52 

0.034 

0.62 

(0.0025) 

0.43 

0.086 

0.53 

0.028 

0.63 

(0.0022) 

0.44 

0.090 

0.54 

0.022 

0.64 

(0.0020) 

0.45 

0.093 

0.55 

0.016 

0.65 

(0.0018) 

0.46 

0.096 

0.56 

0.012 

0.66 

(0.0017) 

0.47 

0.095 

0.57 

0.0079 

0.67 

(0.0016) 

0.48 

0.089 

0.58 

(0.0057) 

0.68 

(0.0015) 

0.49 

0.074 

0.59 

(0.0044) 

0.69 

(0.0014) 

0.70 

(0.0014) 

‘Values  in  parenthesis  are  extrapolated. 
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Fig.  6.  Reflectance  of  {optically  infinite)  deep  water. 

Assuming  the  reflectance  of  the  (optically  infinite)  deep  water  to  be  diffuse*,  the  radiance  reflected  in 
all  upward  directions,  Nf  is 


N  a  =  Hd(0,0)  — 


(2.4.10) 


This  radiance  strikes  the  water-air  surface  and  reflects  at  angle  0.  Some  of  the  radiance  is  transmitted 
through  the  surface,  the  proportion  reflected  is  r(0)  [r(0)  =  r(0)  where  the  relationship  between  6  and  0  is 
defined  by  Snell's  Law,  see  Fig.  7], 


Fig.  7.  Upwelhng  radiance  reflected  at  0  by  amount  rtf?). 


*  Strictly  speaking,  the  reflet  lame  is  not  diffuse  but  m  t  enough  data  on  the  directional  reflectance  is  available  to 
allow  the  reflectance  to  lie  handled  as  other  than  diffuse.  See  Kef  10.  Fig  2(i 
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The  downward  directed  radiance  N(@,<£)  would  thus  be 


N(0,  <f>)  =  Nfr(0) 

The  downwelling  irradiance  from  the  first  water-air  surface  reflectance  would  be: 

HD1(0,0)  =  f2JTN(®,  <j>)  cos®  dtt' 


(2.4.11) 


(2.4.12) 


Substituting  from  Eq.  2.4.10  and  2.4.11. 


HdIO-OR* 

Hm(0,0)  - - J2n  r(0)  cos©  dQ’ 


(2.4.13) 


The  reflectance  of  the  diffuse  upwelling  light  by  the  water-air  interface  is  the  emergent  reflectance  or 
internal  reflectance,  r„ .  Thus, 


i  r 

i",  =  —  /  r(0)  cos©  dfi' 

77  J 


(2.4.14) 


Using  a  discrete  summation  method,  with  Ad  =  5°  and  A<£  =  10°,  rt  has  been  evaluated  to  be  0.485  for  the 
photopic  index  of  refraction  of  sea  water  (n1  =  1.341).  This  is  within  1  percent  of  the  value  reported  by 
Judd12  (n1  =  1.34,  rt  =  0.480)  who  used  a  smaller  summative  interval  (A0  =  2.3°  except  for  the  interval 
near  the  critical  angle).  The  emergent  reflectance  for  fresh  water  and  sail  water  for  all  wavelengths  in  the 
visible  spectrum  varies  less  than  1  percent  from  the  photopic  value  for  sea  water. 

The  downwelling  irradiance  from  the  first  water-air  surface  reflectance  can  now  be  written  thus: 


HD1(0,0)  =  Hd(0,0)  R^r. 


(2.4.15) 


The  downwelling  irradiance  from  the  second  water-air  surface  reflectance  would  be 


HD2(0,0)  =  Hol(0.0)R =  H^O.OHRj, 


(2.4.16) 


Similarly  the  third  and  succeeding  surface  reflectances  would  be 

HD3(0.0)=  Hd(0,  0)  (R00ri)3 


(2.4.17) 


HonfO.O)  =  H^O.OHR^r,)" 


(2.4.18) 


The  radiance  which  is  in  the  direction  which  will  be  refracted  toward  the  sensor's  path  of  sight  6  is: 


N(0p.  A)  -  Hd(0,0)-  +  HD1(0,0)~  +  HD2(0,0)—  +... 


(2.4.19) 


Let  R  r  =  x,  then 

oc  S  ' 


'•S'*®** 


N(0p.0p)  -  Hd<0,0)  — ■  +  Hd(0,  0)  (R«,rm)  +  Hd(0.0)~  (R«rs)2+...  (2.4.20) 


N(0P,  <f>\  =  Hd(0,  0)  —  (’>  +  x  +  x2  +  x3  +  . . .) 

!T 


(2.4.21) 


The  sum  of  the  series  in  the  parenthesis  is  1/(1  -x),  therefore 


N(0p,  <6p)  =  Hd(0,0) 


*d~r.RJ 


(2.4.22) 


Substituting  Eq.  (2.4.5)  in  the  above  equation. 


j  f27rN(0B,  <f>B)  t«?B)  cos0BdflB  +  N,(0,,.O)  t(0.)  cosdt  d:>,  R 
B  ,6)=L‘ - J 


N(0„.  <f>B) 


"O-r.fU 


(2.4.23) 


The  equivalent  derivation  of  the  upwelling  luminance  beneath  the  water  surface  is  as  follows 


Bf  =  ED(0,0)Ro 


(2.4.24) 


Note  the  n  term  is  missing  since  the  luminances  are  in  ft-L. 


8(0.0)  =  Bfr(0) 


(2.4.25) 


1  f 

,(0,0)  =  —  y27rB(0.  0)cos0dfi'  =  Ed(0, 0)  R^r, 


(2.4.26) 


ED2(0,0)  =  Ed(0,  0)  (R^rJ2 


(2.4.27) 


B(0p.0p)  =  Ed(0.  0)  RpjE  Edi(0,  0)  R^  +  ED2(0. 0)  R^.+ . 


(2.4.28) 
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and 


B(6„ 


^ )  = 


Ed{0,  0)  +  Ed(0,  0)  RJfU.)  +  Ed(0,  0)  R00(R00r,)2  +  . . 


(2.4.29) 


and 


B(©p.  *„>  =  Edv0,0) 


R 


00 


O-r.RJ 


(2.4.30) 


Note  the  n  term  is  missing  in  Eq.  2.4.30,  but  when  Eq.  2.4.6  is  substituted,  Eq.  2.4.23  and  2.4.31  become 
equivalent  with  no  n  term  differences. 


B(0P.  <f>P)  =  [/„*..  <£B>  Mb>  c°s^b  dnB  +  B.(0.,O)  t (0.)  cos0.  dft,]Ro 

^H-r,  R„) 


(2.4.31) 


2.4.3  UPWELLING  RADIANCE  (LUMINANCE)  ABOVE  THE  WATER  SURFACE 

The  upwelling  radiance  just  beneath  the  water  surface,  N(0p,<£p),  is  transmitted  through  the  surface, 
and  refracted  by  that  surface  so  that  the  same  energy  is  spread  over  a  larger  solid  angle.  The  projected 
area  of  the  radiance  is  smaller,  being  a  function  of  cos0o  in  the  air,  and  of  cos0o  in  the  water.  The  total 
effect  is  to  reduce  the  radiance  by  the  Fresnel  transmittance  t(0o)  and  by  the  ratio  of  the  respective  solid 
angles  and  projected  areas,  see  Fig.  8. 


Fig.  8.  Upwelling  radiance  at  interface. 
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Nu(0p.<£p)  =  N(©p.0p)  t(0o) 


ds  cos0o  dft' 
ds  cos0o  dQ 


(2.4.32) 


The  solid  angle  of  the  radiance  dfi  may  be  written  in  the  form  sin0od0od<£.  The  azimuth  d0  does  not 
change  above  and  below  the  surface  since  <f>  is  defined  by  the  plane  through  the  light  ray  and  the  normal  to 
the  surface.  Thus: 


Nu(0p,<A  )  =  N(0  ,  <f>  )  t(0o) 


cos0„  sin0„  d0 


cos£?0  sin0o  d0o 


To  obtain  d0  Snell's  law  Eq.  2.1.1  is  differentiated 


0„  =  sin- 


'  7  Sin*°] 


by  proper  substitution. 


d0_  = 


V,_6) sin!e- 


—  cosfldfl 

I  O  O 

n 


n  cos0o 

d@o  =  -  -  d0# 

n  cos0„ 


then  since  s<n0o/sin0o  =  n/n1,  Eq.  2.4.33  becomes 


NU(0P.<£P>  =  N(0p,d.p)  t(0, 


(2.4.33) 


(2.1.1) 


(2.4.34) 


(2.4.35) 


(2.4.36) 


The  change  in  radiance  at  an  interface  due  to  the  square  of  the  ratio  of  the  refractive  indices  was  first 
noted  by  Gershun.13 

The  photometric  expression  for  the  upwelling  luminance  is  directly  comparable  to  Eq.  2.4.36. 


2.5  Directional  Radiance  (Luminance)  of  Calm  Water 

The  directional  radiance  of  calm  water  N((?pfd>p)  is  the  sum  of  the  reflected  sky  and  'or  sunlight  and 
the  upwelling  radiance. 


N,ky  or  ,un(0p.<V  +  Nu(()p,6p) 


Similarly,  for  the  photopic  case 


B «>„.<£„)  =  Blky  or  mun(0p,  6p)  +  Bu(0p,  6p) 


(2.5.1) 


(2.5  2) 
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3.  WIND-RUFFLED  SEA,  AVERAGED 

The  wind-ruffled  sea  has  a  myriad  of  tiny  wave  facets  which  reflect  and  transmit  the  sun  and  sky  light. 
The  time-  or  space-averaged  radiance  (luminance)  of  the  surface  of  the  wind-rumpled  sea  is  dependent 
upon  the  statistics  of  the  tiny  capillary  wave  slopes.  The  statistical  distribution  of  these  wave  facets 
will  be  examined  in  the  first  section.  The  representation  of  the  slope  characteristics  of  the  wind  blown 
sea  by  means  of  an  ergodic  cap  will  be  described  in  the  second  section.  Third,  the  occlusion  from  the  sky 
or  sun  of  a  given  point  on  the  surface  of  the  sea  by  other  waves  is  analyzed.  Then  come  sections  on  the 
equations  for  sun  glitter,  sky  reflectance,  upwelling  light,  and  their  sum  total,  the  time-  or  space-averaged 
directional  radiance  (luminance)  of  the  wind-disturbed  ocean  surface. 

3.1  Wave  Slope  Statistics 

The  distribution  of  the  water  wave  slopes  is  a  function  of  local  wind  speed.  The  relationship  of  the 
slopes  to  the  wind  speed  has  been  measured  by  Duntley14  mostly  at  a  fresh  water  lake  and  by  Cox  and 
Munk15  in  the  open  ocean.  In  both  cases  the  distribution  was  elliptical  with  the  up-wind  direction  con¬ 
taining  larger  slopes  than  the  cross-wind  direction.  The  relationship  of  the  standard  deviation  of  the  wave 
slopes  to  the  wind  speed  is  given  below  for  up-wind  and  cross-wind  for  the  open  ocean.  The  wind  speed, 
v,  is  in  knots  and  is  for  a  41  ft.  height,  the  anemometer  height  normally  used  by  ships  at  sea. 

=  0.00162  f  —  V  (3.1.1) 

\n.mi./ 

a\  =  0.000990  (  ~~~  jv  (3.1.2) 

\n.mi./ 

O 

The  standard  deviation  for  any  azimuth  from  the  wind  direction  Mother  than  up-wind,  4>=  0,  or  cross-wind, 
$  =  90°,  can  be  found  by. 


u\ ■=  - .  (3.1.3) 

*  a 2  sin2  <&+  a\  cos2  <l> 

In  order  to  simplify  the  calculation  of  the  sea  radiance  (luminance)  the  wave  slope  distribution  will  be 
approximated  by  a  circular  rather  than  an  elliptical  probability  function.  This  means  that  the  slope  distri¬ 
bution  in  the  up-wind  and  cross-wind  directions  will  be  assumed  to  be  equivalent.  In  doing  so,  to 
minimize  the  error,  the  standard  deviation  for  a  45°  angle  from  the  wind  will  be  used  since  this  lies  be¬ 
tween  the  measured  values  for  up-wind  and  cross-wind. 


Also, 


a*  +  ai 


0.00123 


(3.1.4) 


a  -  0.0351  V7. 


(3.1.5) 
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This  means  that  the  wave  slope  distribution  will  be  represented  three  dimensionally  as  illustrated  in 
Fig.  9. 


y 


Where  the  height,  y  :s 


and  the  radius  X  is 


X  = 


tan# 


a 


(3.1.6) 


(3.1.7) 


tan 0  is  the  wave  slope.  The  volume  V  of  the  solid  in  Fig.  9  is  the  probability  of  the  wave  slopes.  Thus, 
the  total  volume  or  total  probability  V  =  1.  From  Fig.  9. 


AV  =  yir(X*  -  Xj)  =  iryiX2  -  X,)  (X2  +  X,) 


(3.1.8) 


and 


A  V  =  2  ir  y  X  A  X 


tl-  i 


(3.1.9) 


(3  1.10) 


15 


3.2  Ergodic  Cap 

The  slope  distribution  of  the  wave  facets  for  a  given  wind  speed  can  be  represented  by  a  three- 
dimensional  ergodic  cap.  The  surface  of  the  cap  is  shaped  so  that  it  has  the  correct  distribution  of  slopes 
for  that  wind  speed.  The  ergodic  cap  is  circular  since  the  probability  function  has  been  assumed  to  be 
circular.  The  horizontal  projected  area  AH  of  the  ergodic  cap  is  equated  to  the  probability  or  volume  un¬ 
der  the  three-dimensional  probability  function  illustrated  in  Fig.  9.  The  total  horizontal  projected  area 
thus  equals  1. 

*The  derivation  for  a  non-circular  ergodic  cop  bus  been  made  by  Preisendorfor  in  Kef.  2.  The  cap  becomes  elliptical 
instead  of  circular  but  the  general  form  of  the  equations  is  .similar. 
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- IV  -  '  •  '***» 


Ah  =  ;rR2  =  1  (3.2.1) 

The  portion  of  the  ergodic  cap  which  would  have  slopes  between  tan0,  and  tan02  in  all  azimuths 
would  be  between  r,  and  r2,  see  Fig.  10.  where, 

(3.2.2) 

(3.2.3) 


Ah(A  0)  =  V(A0) 


2  2  (  tan202\  (  tan20i 

77 (r2  -  r.  )  =  -  exp  ( - x  J  +  exp  I - 

\  2°  J  \  2 a2  , 


Also,  the  portion  of  the  ergcdic  cap  with  slopes  between  tan0,  and  tan02  and  with  normals  from  the  sur¬ 
face  lying  between  <f>l  and  <62  would  be, 


AH(A0,A<3)  =  V(A0,A<3) 


(3.2.4) 


A  4> 


/  2  2  > 

<f2  "  r,  >  =  — 

2  77 


tan2fO 


-  exp  - 


2a2 


+  exp 


tan 2  0, 


2a2 


(3.2.5) 


The  true  area  of  the  surface  of  the  cap  with  normals  from  the  surface  lying  between  0,  and  02,  and  <3, 
and  <£2.  is 


Ah(A0,A  6) 

At(A0,A<3)=  - 

cos0N 


where  0N  is  (0,  +  02)/ 2  or  the  average  slope  of  the  elemental  surface  area. 


(3.2.6) 
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The  area  as  seen  by  an  observer  or  sensor  along  an  inclined  path  of  sight  at  0p,  tf>  ,  is  equal  to  the 
true  area  projected  onto  the  plane  perpendicular  to  the  path  of  sight.  This  projected  area  is  equal  to  the 
true  area  times  the  cosine  of  the  angle  between  the  normal  from  the  surface  at  0N,  <£N  [where  <pN  = 
|(£j  +  <^)/2]  and  the  observer's  position,  at  0o,  <f>o.  This  angle  is  the  Fresnel  reflectance  angle,  0R. 

\{A0,A<f>)  =  AT(Atf,A0)  cos0R  (3.2.7) 


When  0,  r  go0,  Ap(A6,A(f>)  =  0. 

The  elemental  horizontal  area  of  the  ergodic  cap,  AH(A0,  A<£),  is  equivalent  to  the  probability  of  oc¬ 
currence  of  the  wave  slopes  between  AO  a.id  A <f>.  The  elemental  projected  area  of  the  ergodic  cap  divided 
by  the  total  projected  area  of  the  ergodic  cap,  A5(A0,  A<t>)/Ap,  is  the  probability  of  occurrence  of  the  radi¬ 
ance  (luminance)  reflected  and  transmitted  by  the  elemental  surface  area  in  the  direction  of  the  path  of 
sight  of  the  observer. 

The  erqodic  cap  will  be  used  to  develop  the  geometry  of  the  reflectance  and  transmittance  of  the  sun, 
sky,  and  upwelling  light.  The  probability  function  will  be  used  to  give  the  probability  of  occurrence  of  the 
wave  slopes.  The  ergodic  cap  will  be  used  as  a  continuously  curved  surface  for  development  of  the  equa¬ 
tions  for  reflectance  and  transmittance  of  the  sun  by  the  water  surface.  The  reflectance  and  transmittance 
of  the  sky  light,  however,  will  be  handled  by  a  simplified  ergodic  cap,  having,  instead  of  a  continuously 
curved  surface,  a  series  of  flat  surfaces,  each  having  a  normal  from  the  surface  at  angle  0N,<£N. 

3.2.1  ROTATION  OF  COORDINATE  SYSTEMS 

Computation  of  the  reflectance  (transmittance)  angle  for  each  ergodic  can  element,,  and  the  position 
0B,  <£Bof  the  sky  radiance  (luminance)  reflected  (oi  transmitted)  is  facilitated  by  using  unit  vector  nota¬ 
tion,  the  normal  from  the  surface  being  N,  the  sky  position  being  B,  and  the  observer  position,  0.  By  ro¬ 
tating  the  coordinate  system  from  the  zenith,  sun  azimuth  system  to  the  coordinate  system  defined  by  the 
observer  position,  the  reflectance  angles  and  sky  position  can  be  found  easily.  Coordinate  rotc.ion  is 
also  useful  for  finding  the  portion  of  the  ergodic  cap  which  reflects  the  sun  light  into  the  path  of  sight  of 
the  observer. 

The  derivation  of  the  equations  for  the  transformation  of  rotation  nf  axis  expressed  in  spherical  co¬ 
ordinates  is  given  as  Appendix  A.  The  angles  in  the  unrotated  syste..,  are  given  as  unpiimed:  zenith 
angle  0.  azimuth  from  sun,  <f>.  See  Fig.  11. 

k 


Fig.  11.  Unrotated  coordinate  system. 
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The  rotated  coordinate  system  is  defined  by  angles  0t  and  <£,  in  the  unrotated  system.  See  Fig.  12. 


«  k' 


Fig.  12.  Rotated  axes. 


The  angles  in  the  rotated  system  are  denoted  as  primed  values,  O'  <£'•  See  Fig.  13. 


The  baL.c  equations  for  finding  direction  angles  in  a  rotated  system  are  as  follows  (equations  from 
Group  A.  12  in  Appendix  A). 

cost?'  =  sinflcos^sinf^  cos^  +  sinflsin^sinf?,  sin^,  +  cosdJcosflj  (3.2.8) 

cos<£'  =  tsin^cos^sin^cos^),  +  sin0sin<£cos01sin<£1  -cos^sin^]  *  sin0‘  (3.2.9) 


and 


sin<£'  =  [— sin0cos0  sin^,  +  sin^sin^cos^]  -  sin0'  (3.2.10) 

Eqs.  3.2.8  through  3.2.10  are  for  finding  the  zenith  angle  O'  and  azimuth  <£'  in  the  rotated  system  for  the 
vector  at  zenith  angle  6  and  azimuth  <f>  in  the  unrotated  system.  The  rotated  axes  are  defined  by  zenith 
angle  0X  and  azimuth  <£r  It  is  necessary  to  find  both  sinks'  and  cos<£'  in  order  to  place  <f>'  in  the  correct 
quadrant. 

The  basic  equations  for  finding  direction  angles  in  the  unrotated  system  are  as  follows  (equations  are 
from  Group  A.11  in  Appendix  A). 


cos0  =  -  sin^'cos^' sin0j  +  cos0’ cos0,  (3.2.11) 

cos <f>  -  [sin0'  cos^1  cos0t  cos^j  -sin0'  sin<£'  sin^  +  cos0'  sin0j  cos<£,]  t  sin0  (3.2.12) 

sin<£  =  [sin0‘cos^'cos01  sin^,  +  sin0‘ sin^'cos^  +  cos0'  sin0,  sin<£,]  -  sin0  (3.2.13) 

Eqs.  3.2.11  through  3.2.13  are  for  computing  the  zr.iith  angle  0  and  azimuth  <f>  in  the  unrotated  system  for 
the  vector  at  zenith  angle  0'  and  azimuth  cf>'  in  the  rotated  system.  The  rotated  axes  are  again  defined  by 
zenith  angle  0,  and  azimuth  <f>l.  Again  i<  is  necessary  to  find  both  sin$  and  cos<f>  in  order  to  place  the 
azimuth  <f>  in  the  correct  quadrant. 


3.3  Occupation  from  Light  Sources 

The  ergodic  cap  simulates  the  slope  characteristics  of  the  wind-ruffled  sea,  but  does  not  simulate  the 
effect  of  hiding  of  one  wave  by  another.  A  given  point  on  the  surface  of  the  sea  will  sometimes  be  ob¬ 
scured  from  a  light  source;  this  effect  increases  as  the  zenith  angle  of  the  light  source  becomes  large. 
Q°  (*10  is  the  average  fractional  time  that  an  incident  source  from  elevation  angle  breaches  a  given  point 
on  the  sea  surface  without  f>st  having  intersected  the  water  surface  at  some  other  point.  The  derivation 
of  Q°(*P)  is  given  in  Appendix  B  by  MacAdam,  it  is  consistent  with  the  derivations  given  in  Reference  2. 


a°m  =  exp - -f— 

tan  'VyfT-n  _ 


(3.3.1) 


Changing  from  elevation  angle  notation  to  zer:»h  angle  notation,  and  substituting  from  Eq.  3.1.5  for  a, 
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Q°(0)  -•  exp  (  -  0.0140  tan 0  /7h 


(3.3.2) 


Q°  (0)  is  the  fraction  of  the  light  from  a  source  at  zenith  angle  0  which  reaches  the  sea  surface. 


3.4  Sun  Glitter 

The  time-averaged  sunlight  reflected  by  the  water-wave  facets  is  called  the  sun  glitter.  The  average 
apparent  sun  radiance,  Ns  (0s,  0)  is  partially  hidden  from  the  sea  surface  by  other  waves.  The  radiance 
that  reaches  the  water  is  Ns  (0s,  0)  Q °(0s).  This  radiance  is  reflected  at  angle  0SR  into  the  path  of  sight 
of  the  sensor.  The  proportion  of  the  time  this  sun  radiance  will  be  reflected  into  the  path  of  sight  is 
Ap(At)/Ap  Ap  is  the  total  projected  area  of  the  ergodic  cap.  The  equation  for  Ap  will  be  derived  in  the 
next  section  which  will  concern  the  sky  reflectance  for  the  wind-ruffled  sea 

Ng(0p.<£p)  =  NS(0SO)  Qo(0a)  r  (0SR)  Ap(At)/Ap  (3  4.1) 

The  photometric  equation  is  directly  analogous. 


Bg(0 P,6P)  =  8a(0s,O)  Q°(0a)  r  (0SR)  Ap(At)/Ap 


(3.4.2) 


Ap(At)  is  the  portion  of  the  projected  ergodic  cap  that  is  tilted  at  the  appropriate  angles  to  reflect  the  sun 
into  the  path  of  sight.  Thus,, 


Ap<  At) 


V(At)  cos0SR 
cos0SN 


(3.4.3) 


where  V(At)  is  the  probability  of  the  wave  slopes  which  reflect  the  sun  into  the  path  of  sight,  0gR  is  the 

angle  between  the  observer  and  the  normal  from  the  facet  reflecting  the  center  of  the  sun,  and  0  is  the 

zenith  angle  of  the  normal  from  the  facet  reflecting  the  center  of  the  sun,  the  probability  V(At)  will  be  de¬ 

rived  first,  and  then  the  equations  for  computing  the  angles  0SR  and  0SN  will  be  set  forth. 

The  first  step  in  deriving  the  probability  will  be  to  redraw  the  probability  function  from  Fig.  9  by 
changing  variables,  see  Fig.  14.  The  radius  changes  from  X  to  tan0  and  the  height  y  becomes  y'  and  is 

1  /  tan20  \ 

y'  =  —  exp  ( - 1  (3  4.4) 

2 no1  \  2n2  / 
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V' 


1  tea  (  ta"20\ 

V  = —  Jo  tan0expl-__  ldtan0  (3.4.5) 


=  1 


(3.4.6) 


The  element  of  volume  of  the  probability  function  representing  the  wave  slopes  which  reflect  the  solar 
disk  of  size  ft.,  can  reasonably  be  represented  as  an  elliptical  cylinder.  The  center  of  the  cylinder  is  at 
tan0SN,  the  slope  which  reflects  the  center  of  the  sun.  The  height  of  the  cylinder  is  y'a  where 


1 


tan20sfA 

2c72  )  ■ 


(3.4.7) 


The  ellipse  of  the  cylinder  is  the  tolerance  ellipse  as  defined  by  Cox  and  Munk  Ref.  14.  Its  area  At 
is  given  as 


At  =  'An  t2  sec3  /3o  sec«o 


(3.4.8) 


where  e  is  the  angular  radius  of  the  sun  in  radians,  /3o  is  the  tilt  angle  and  is  equivalent  to  f"*  angle  0SN 
in  the  notation  herein,  and  <»0  is  the  angle  of  reflection  used  herein  as  0SR.  Thus  Eq.  3  4.8  can  be  rewritten, 
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(3.4.9) 


At  =  'A  rr  ( 2  sec30SN  sec0SR. 

The  angular  radius  of  the  sun  for  mean  solar  distance  is 

e  =  0.004659  radians  (3.4.10) 

The  volume  of  an  elliptical  cylinder  is  the  height  times  the  area  of  the  ellipse.  Thus  combining  Eqs.  3.4.7 
and  3.4.9, 


V(At)  = 


<-2sec30SN  sec0SR 
Ba2 


tan  2°sn  \ 
2a2  )  ' 


(3.4.11) 


Finally  combining  Eqs.  3.4.3  and  3.4.11,,  the  projected  area  of  the  ergodic  cap  which  reflects  the  sun  into 
the  path  of  sight  is. 


Ap(At)  =  - 

8o2cos40S(J 


Substituting  for  a 2  from  Eq.  3.1.4, 


Ap(At)  - - 

0.00984v  cos40sn 


(3.4.12) 


(3.4.13) 


The  angles  0SR  and  0SN  are  easily  computed  by  rotating  from  one  coordinate  system  to  another  and 
back  again.  The  principal  or  unrotated  coordinate  system  is  defined  by  the  true  zenith,  and  the  azimuth 
from  the  sun.  The  prime  coordinate  system  is  defined  by  a  vector  S  toward  the  sun  at  d3  and  <£s  =  0. 
See  Fig.  15.  k 


23 


The  observers  position  vector  has  a  zenith  angle  0o  and  azimuth  <jSo.  The  relationship  of  the  observer  an¬ 
gles  Qo  ,  6a  to  the  path  of  sight  angles  was  discussed  in  section  2.3.  The  observer’s  position  in  the 
rotated  system  is  found  by  appropriate  substitution  into  Eqs.  3.2  8  and  3.2.9.  Since  <£,  =  <f>s  =  0,  some  of 
some  of  the  terms  drop  out. 


cos^  =  sin$0  cos<£0  sin0s  +  cos0ocos0s  (3.4.14) 

cos<^  =  [sin0ocos<£ocos0s  -cosf?0 sin#,]  *  sin6£  .  (3.4.15) 

It  is  not  necessary  in  this  case  to  determine  directly  only  cos<^ ,  the  reason  for  this  will  be  seen 
shortly.  Therefore,  it  is  not  necessary  to  compute  sin^  from  Eq.  3.2.10.  Since  6 £  is  the  angle  between 
the  sun  S  and  the  observer  0 ,  the  Fresnel  reflectance  angle,  0SR,  is  equal  to  half  of  6£. 


e 


SR  _ 


2 


(3.4.16) 


The  normal  from  the  facet  reflecting  the  center  of  the  sun  SN  is  at  6£/2,  <f>0  in  the  prime  coordinate  system, 
see  Fig.  16. 


Rotating  back  to  the  principal  coordinate  system  by  substituting  appropriately  in  Eq.  3.2.11,  Eq.  3.4.17  is 
obtained. 


el 


e' 


cos0SN  =  -sin  —  cosc^  sin0s  +  cos  —  cos0, . 


(3.4.17) 


3.5  Sky  Reflectance 

The  time-averaged  wind-ruffled  sea  reflects  light  from  various  portions  of  the  sky.  To  simplify  the  cal¬ 
culation.  the  curved  surface  of  the  ergodic  cap  is  represented  by  a  series  of  flat  surfaces.  Thus,  the  time- 
averaged  inherent  radiance  of  the  sea  due  to  sky  reflectance  can  be  calculated  by  summing  the  appropri¬ 
ately  weighted  reflected  sky  radiances  by  each  of  the  flat  surfaces. 
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1 


^sky^P'^p)  =  ~ 


n  f  m 

-*■  -s. 


N(^b'  ‘Wij  r  (^BRij)  C0S^BRIj 


J 


Ah(«m.) 


COS0K 


(3.5.1) 


The  equation  is  similar  for  the  photopic  case. 

The  surface  of  the  ergodic  cap  can  be  cut  into  sections  represented  by  flat  surfaces.  Each  flat  surface 


is  represented  by  a  normal  from  the  surface  at  dN  and  where  6^  =  (0,  +  02)/ 2  and  <f>u  =  (^  +  <f>2)/ 2; 

A0n  =  02  -  0,  and  A<£N  =  <j>2  -  <f>t.  The  slopes  from  0N  =  0°to  A0/2  regardless  of  azimuth  are  repre¬ 
sented  by  a  horizontal  surface,  0N  =  0.  Thus,  from  Eqs.  3.2.2  and  3.2.3 


/  tan2(A0K/2)' 

Ah(0)  =  -  exp  I - : - I  +  1 


2(7  2 


(3.5.2) 


and  from  Eqs.  3.2.4  and  3.2.5 


1  f  /  tan2(0N  +  A0N/2)  \  /  tan2(0N  -  A0N/2)  \' 

rxpv — — ;+en- — & — ) 


(3.5.3) 


The  term  1/m  =  A<f>/2n  where  A 0  is  in  radians.  The  number  n  of  i  or  0N  terms  depends  upon  wind  speed 
and  the  assigned  precision  level  for  the  total  wave  slope  probability.  The  number  of  0N  terms  for  an  er¬ 
godic  cap  with  A0N  =  5°,  with  slopes  greater  than  0Nn  +5° contributing  less  than  10"4  to  the  total  wave 
slope  probability  is  given  in  Table  3  for  various  wind  speeds. 


Table  3.  Ergodic  cap  specifications  for  A  0N  =  5°. 


Wind  speed  (knots) 
8  14  20 


°Nn  (degrees) 


3 

10 


4 

15 


6 

25 


7 

30 


8 

35 


of  the  ergodic  rap 


AP(0N)  = 


cos0BRAH(0N) 


COS0* 


Thus,  the  total  projected  area  is  the  sum  of  all  these  elemental  projected  areas 

A„tfU 


n  ~  m  ~  'N  ' 

-  cos0BR.  - 

1  ~  1  >  1  cos 0 

L.  N  i 


•s 


I 


& 

■'*£ 


26 

32 

38 

9 

9 

10 

40 

40 

45 

of  each 

flat  surface  or  element 

(3.5.4) 

areas. 

(3.5  5) 
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The  total  projected  area  is  within  1  percent  of  the  cosine  of  the  observer  angle,  Ap  =  cos 0o,  for  observer 
zenith  angles  0  to  70°.  For  larger  observer  zenith  angles  the  projected  area  becomes  a  function  of  wind 
speed.  Evaluation  of  the  projected  area  for  large  observer  zenith  angles  for  an  ergodic  cap  with  flat  sur¬ 
faces  representing  A0N  =5°  and  A d>N  =  10°,  is  given  in  Table  4. 


Table  4.  Total  projected  area,  Ap,  of  ergodic  cap  for  large  observer  zenith  angles. 


Observer 

7pnith 

Wind  Speed  (Knots) 

i-CI  II  U  1 

Angle  0O 

0 

2 

4 

8 

14 

20 

26 

32 

38 

75 

0.259 

0.259 

0.259 

0.259 

0.260 

0.262 

0.264 

0.267 

0.270 

80 

0.174 

0.174 

0.174 

0.175 

0.179 

0.184 

0.189 

0.194 

0.199 

85 

0.0872 

0.0878 

0.0907 

0.0977 

0.107 

0.116 

0.124 

0.131 

0.137 

90 

0.00000 

0.0198 

0.0282 

0.0398 

0.0527 

0.0630 

0.0718 

0.0796 

0.0868 

The  reflectance  angle  0BR  and  the  sky  position  angles  0B,  cf>B,  can  be  computed  by  assuming  unit  vec¬ 
tor  notation  and  rotating  axes.  The  observer  position  vector  0  at  6o,  <j>o  defines  the  rotated  system.  The 
normal  N  from  a  given  ergodic  cap  element  is  at  0N,  <£N.  The  first  task  is  to  find  the  normal  direction  an¬ 
gles,  0N\  <£n\  in  the  rotated  system.  Appropriately  substituting  into  Eqs.  3.2.8,  3.2.9,  and  3.2.10. 


cos0R  =  sin0Ncos<£Nsin0ocos0o  +  sin0Nsin0Nsin0osin<j!.o  +  cos0Ncos0o  (3.5.6) 

cos<£R  =  [sin0Ncos<£Ncos0ocos<£o  +  sin0Nsin<£Ncos0osin0o  -  cos0Nsin0j  *  sin0R  (3.5.7) 

sin^  =  [- sin0Ncos<^>Nsin9So  +  sin0Nsin<£Ncos0j  *  sin0R  (3.5.8) 


0^  is  also  the  Fresnel  reflectance  angle,  0BR.  When  0^  =  0BR>  90°,  the  reflectance  is  zero  and  the  pro¬ 
jected  area  is  zero.  The  sky  position  in  the  rotated  system  is  twice  the  Fresnel  reflectance  angle  and  in 
the  same  azimuth  as  the  normal  from  the  reflecting  surface.  Thus 


0B  =  20^  (3.5.9) 

4>B  =  (3-5.10) 


The  sky  position  in  the  unrotated  system  is  determined  by  appropriately  substituting  into  Eqs.  3.2.11 
through  3.2.13.  See  Fig.  17. 
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cos0B  =  sin0Bcos0Bsin0o  +  cos0'^cos0o  (3.5.11) 

cos 4>b  ~  [sin0Bcos<£Bcos0ocos0o  -  sin0Bsin0gSin(£o  +  cos0Bsin0ocos<£o]  *  sin0B  (3.5.12) 

sin^B  =  [sin0Bcos<£Bcos0osin<£o  +  sinfl^in^gCos^  +  cos0Bsin0osin<£j  -  sin0B  (3.5.13) 

When  0b>9O°,  the  sky  radiance  (luminance)  is  considered  =  zero  since  the  contribution  from  the  lower  sky 
(sea  surface)  is  negligible.  Also,  when  the  sky  luminance  distribution  is  symmetrical  with  respect  to  the 
azimuth  of  the  sun  (clear  days  without  clouds  are  reasonably  symmetrical)  it  is  not  necessary  tc  compute 
sin^ig;  cos0B  will  yield  answers  between  0  and  180°  in  azimuth. 


3.6  Upwelling  Light 


The  light  upwelling  from  the  surface  of  the  wind-ruffled  sea  is  similar  in  form  to  the  light  upwelling 
from  a  calm  sea,  (see  Eqs.  2.4.23  and  2.4.36),  with  these  differences:  Each  source  radiance,  N(0B,  $B), 
will  be  accompanied  by  an  occultation  term,  Q°(6B),  and  a  time-averaged  transmittance  and  cosine, 
/t(0B)cos0Bdt  for  the  transmittance  and  redirection  of  the  light  ray  beneath  the  surface  of  the  sea  by  the 
tilted  wave  facets.  The  internal  reflectance  of  the  water  r,  will  be  a  time-averaged  internal  reflectance 
of  diffuse  light  /r^dt  Finally,  the  transmittance  of  the  diffuse  light  upward  in  the  direction  of  the  path 
of  sight  is  a  time-averaged  transmittance  /t(0)dt.  Thus, 


N J0P.  <f>J 


=  ^  f2nN(diy  *B>  QO(0Byt(0g)COS0gdtdfiB 
+  N„(0.,  0)Q °{0,)/t(0,)  cos0.dtdQ.j  x 


- 'r~  /^)d?(— j 


(3.6.1) 
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The  photometric  equation  will  not  be  wiitten  out  since  it  is  necessary  only  to  substitute  the  luminance 
symbol,  B,  for  every  radiance  symbol,  N. 


3.6.1  TIME-AVERAGED  TRANSMITTANCE  AND  COSINE 

The  slope  distribution  of  the  wind-ruffled  sea  surface  is  represented  by  an  ergodic  cap  with  flat  ele¬ 
ments  as  was  described  in  the  preceding  section  on  sky  reflectance.  For  a  single  source  at  angle  0B,  the 
occultation  is  invariant  but  the  Fresnel  transmittance  angle,  0gR,  is  different  for  each  flat  element  of  the 
ergodic  cap.  Second,  the  ratio  of  the  solid  angles  and  projected  areas  changes  with  each  element, 
(ds  cos0BRdil)  *  (ds  cos0BRdfl')  (see  Eq.  2.4.1],  Third,  a  flat  plate  collector  will  receive  the  radiance 
at  the  refracted  angle  0B  and  solid  angle  dQ'  which  will  be  different  for  each  element  of  the  ergodic  cap. 
Fourth,  the  probability  of  transmitting  the  radiance  from  B  by  the  elemental  surface  area  represented  by 
the  ergodic  cap  is  cos0BRAH(0N)/Apcos0N.  Thus,  for  one  element  of  the  ergodic  cap,  the  irradiance  be¬ 
neath  the  surface  of  the  water,  AH  (0,0),  received  from  the  source  radiance  at  0B,  4>B  is 

1  r  ds  cosflBRdft  Ah«?n)  n 

AH(0, 0)  =—  N(0B,<£B)Q°(0B)t(0BR)  -  cos0Bdfl'  cos 0- -  I  (3.6.2) 

AP  L  ds  cos0BRdfi'  cos0N  J 


The  irradiance  received  from  the  source  radiance  through  all  elements  of  the  ergodic  cap  is,  as  follows, 


H(0,0)  =  N(0  ,  <j>J  G°(0J  < 


n  /  m  COs2^BRij 

2  (  2  t(0BRij)  - - -  COS0 

‘-‘V"1  '  C0S®BRii  1 


COSdN. 


n  /  m 

2  2 

t.  t \i«  s 


COS0 


BRij 


a„(0ni) 


COS0M. 

W  l 


dfl. 


(3.6.3) 


The  term  in  braces  is  called  the  time-averaged  transmittance  and  cosin e,  J t(0)cos0dt. 


yt(0B)cos0Bdt  = 


2 

i=  i 


cos  20 


^BRij^ 


BRij 


j~  l 


COS0 


COS© 


BiJ 


BRij 


AH(0Ni) 


COS0 


Ni 


n  m 

2  2 

i=  1  j=  i 


COS0 


BRij 


AH(0Ni) 


(3.6.4) 


cos0 


Ni 
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Thus,  the  irradiance  received  from  the  source  radiance  at  is 

B 


H(0,0)  =  N(0. 


<t>B)  Q°(«b> 


]h 


,)  cos0Bdt  dftB 


(3.6.5) 


and  the  irradiance  from  the  upper  hemisphere  is 


h 

J  2J7 


H(0,0)  =  /  n(&b,  0b)  q°(^b) J ne 


cos0Bdt  d{2B 


(3.6.6) 


The  Fresnel  transmittance^angle  0BR  is  the  angle  between  a  vector  in  the  source  diiection,  IT  and  the 
normal  from  the  fiat  surface,  N.  The  axes  are  rotated  to  the  vector  N  and  the  luminance  zenith  angle,  dB, 
determined  in  the  prime  system.  d'B  =  dBR.  See  Fig.  18. 
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Since  the  time-averaged  transmittance  and  cosine  should  not  change  with  source  azimuth,  the  integral  in 
Eq.  3.6.4  can  be  evaluated  by  assuming  ^  =  0.  Thus,  from  Eq.  3.2.8, 


cos^br=  sin(?Br,in0N  coscSN+  cos0Bcos0N 


(3.6.7) 


When  00R  v  go°  both  the  transmittancp  and  ihe  projected  area  of  the  element  are  zero  and  that  element  of 
tne  ergodic  cap  is  not  included  in  the  summation.  The  refracted  angle  0BR  is  found  by  means  of  Snell’s 
law,  Eq.  2.1.1, 


The  zenith  angle  of  the  refracted  ray  ©B  can  be  found  by  noting  the  analogy  between  the  angles  be¬ 
neath  the  surface  of  the  water  and  above  the  sutface  of  the  water.  A  new  vector  §  is  defined  in  the  op¬ 
posite  direction  to  the  refracted  ray,  see  Fig.  19. 


Fig.  19.  Plane  passed  through  N  and  B. 


<f>'B  is  the  azimuth  of  both  the  radiance  vector  B  and  the  refracted  vector  IB  in  the  prime  system  defined  by 
the  normal  N.  The  zenith  angle  of  the  refracted  vector  in  the  prime  system  0B  =  0Bp.  cos0B  is  found  by 
appropriate  substitution  into  Eq.  3.2.9. 

cos<£b  =  £  sin0B  cos0N  cos 0N  -  cos0B  sin0N  J  *  sin(?BR  (3.6.9) 


Now  the  zenith  angle  of  the  refracted  ray  0B  can  be  found  by  rotating  from  the  prime  system  to  the  princi¬ 
pal  system  by  using  Eq.  3.2.11  appropriately. 


cosQB  =  -sin0'B  cos0B  sin0N  +  cos0B  cos0N  (3.6.10) 

The  time-averaged  transmittance  and  cosine  was  evaluated  for  various  source  radiance  angles  and 
wind  speeds  for  an  index  of  refraction  n'  =  1.341  (photopic,  sea  water)  with  an  ergodic  cap  element  size 
of  A 6  =  5°  and  A<£  =  10°.  The  result  is  given  in  Table  5.  These  values  may  also  be  used  for  other  wave¬ 
lengths  in  the  visible  spectrum  for  both  salt  water  or  fiesh  water  if  the  accuracy  desired  is  no  greater  than 
1  percent. 


r 


Table  5.  Time-averaged  transmittance  and  cosine  for  source  at  zenith  angle,  theta. 


Source 

Zenith 

Angle 

0 

2 

4 

Wind 

8 

0 

0.979 

0.977 

0.976 

0.974 

5 

0.975 

0.974 

0.972 

0.970 

10 

0.964 

0.962 

0.961 

0.959 

15 

0.945 

0.944 

0.943 

0.940 

20 

0.920 

0.918 

0.917 

0.914 

25 

0.887 

0.885 

0.884 

0.881 

30 

0.847 

0.845 

0.844 

0.841 

35 

0.800 

0.798 

0.797 

0.794 

40 

0.747 

0.745 

0.743 

0.741 

45 

0.687 

0.685 

0.684 

0.681 

50 

0.620 

0.619 

0.617 

0.615 

55 

0.548 

0.546 

0.545 

0.543 

60 

0.469 

0.468 

0.468 

0.466 

65 

0.385 

0.385 

0.385 

0.386 

70 

0.295 

0.298 

0.299 

0.303 

75 

0.203 

0.209 

0.214 

0.224 

80 

0.113 

0.126 

0.136 

0.153 

85 

0.0361 

0.0610 

0.0751 

0.0969 

90 

0 

0.0265 

0.0390 

0.0594 

Speed  (Knots) 


14 

20 

26 

32 

38 

0.970 

0.967 

0.963 

0.960 

0.956 

0.966 

0.963 

0.959 

0.956 

0.952 

0.955 

0.951 

0.948 

0.944 

0.941 

0.936 

0.933 

0.929 

0.926 

0.922 

0.910 

0.907 

0.903 

0.899 

0.896 

0.877 

0.873 

0.870 

0.866 

0.863 

0.837 

0.833 

0.829 

0.826 

0.822 

0.790 

0.786 

0.782 

0.779 

0.775 

0.736 

0.733 

0.729 

0.725 

0.722 

0.677 

0.673 

0.669 

0.666 

0.663 

0.611 

0.608 

0.605 

0.602 

0.599 

0.540 

0.538 

0.536 

0.534 

0.532 

0.465 

0.464 

0.464 

0.464 

0.463 

0.387 

0.389 

0.391 

0.393 

0.395 

0.310 

0.315 

0.321 

0.325 

0.329 

0.236 

0.247 

0.255 

0.262 

0.268 

0.172 

0.186 

0.197 

0.206 

0.213 

0.119 

0.135 

0.148 

0.157 

0.165 

0.0809 

0.0961 

0.108 

0.117 

0.124 

3.6.2  TIME-AVERAGED  INTERNAL  REFLECTANCE  OF  DIFFUSE  LIGHT 

For  a  wind-ruffled  sea,  the  equation  for  the  internal  reflectance  of  diffuse  light,  Eq.  2.4.14,  changes  in 
the  following  manner 


r  -  1  r 

C 

Kelt  =- 

1  r(0)  dt 

J  rr  J  2Tr 

J  j 

cos©  dO1. 


(3.6.11) 


The  term  in  the  brackets  is  the  time-averaged  internal  reflectance  of  diffuse  light  into  the  path  of  sight 
0  of  a  sensor.  See  Fig.  20. 


For  this  computation  the  ergodic  cap  is  assumed  to  have  a  normal  N'  from  beneath  the  surface  at  a  nadir 
angle  ©N,  and  azimuth,  <1^  where  in  this  instance  0N  =  0N  and  =  <£N  +  180°.  See  Fig.  21. 


Fig.  21.  Element  of  ergodic  cap. 

Since  the  time-averaged  ref'ectance  should  not  change  with  sensor  azimuth,  $o  is  assumed  equal  to  zero. 
The  reflectance  angle  ^  is  found  by  appropriate  substitution  into  Eq.  3.2.8. 

cos©R  =  sin0N  cos<J)N  sin0o  +  cos0N  cos0o  (3.6.13) 


The  angle  0R  is  found  from  Snell's  law,  rephrasing  Eq.  2.1.1, 


(3.6.14) 


32 


^V*-**"**®*?'  *#*'■**« 


” ”  >'  ’ w  wip^  /  • 11  r‘ 


•./SWt**'*'' 


The  Fresnel  reflectance  r(0R)  is  found  directly  from  Eq.  2.1.2  r(0R)  =  r(0R).  The  total  projected  area  of  the 
ergodic  cap  is  the  same  as  before,  or  rephrased  in  terms  of  nadir  angles 


2  2  COS0, 


1=11=1 


AH(e„.) 

cos0Ni 


(3.6.15) 


The  time-averaged  internal  reflectance  for  diffuse  light  was  evaluated  for  an  index  of  refraction 
n‘  =  1.341  (photopic  sea  water)  with  an  ergodic  cap  element  size  of  A0  =  5°,  A$  -  10°.  The  summative 
grid  for  Eq.  3.6.12  was  the  same.  The  result  is  given  in  Table  6  for  wind  speeds  0  to  38  knots.  These 
values  may  be  used  for  fresh  water  as  well  as  salt  water,  for  all  wavelengths  in  the  visible  spectrum  with 
an  accuracy  or  2  percent. 


Table  6.  Time-averaged  internal  reflectance  of  diffuse  light  for  n'  =  1.341 


Wind  Speed  (Knots) 


20 

26 

32 

38 

0.470 

0.466 

0.463 

0.459 

3.6.3  TIME-AVERAGED  TRANSMITTANCE 


The  time-averaged  Fresnel  transmittance  toward  0o  of  a  uniform  radiance  (luminance)  from  beneath  the 
surface  of  the  water  is  equivalent  to  the  time-averaged  transmittance  from  air-to-water  from  a  source  at  0o. 
For  ease  in  computation  let  0B  =  0o.  Then,  the  time-averaged  transmittance  is  the  sum  of  the  transmit- 
tances  of  each  of  the  elements  of  the  ergodic  cap  weighted  by  the  probability  appropriate  to  the  source 


direction  0B. 


t(0R)  dt  =  —{22  t(0nRii)  cost? 


[  i=  i  i=i 


cosf?„ . 


(3.6.161 


A  is  the  projected  area  as  defined  in  Eq.  3.5.6.  Since  the  time-averaged  transmittance  should  not  change 
with  source  azimuth,  the  integral  is  evaluated  by  assuming  =  0.  Thus  0BR  can  be  found  from 
Eq.  3.6.7. 
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The  time-averaged  transmittance  has  been  evaluated  for  the  index  of  refraction  n'  =  1.341  (photopic  in¬ 
dex  for  sea  water).  The  ergodic  cap  was  again  made  up  of  flat  surfaces  representing  A(?N  =  5°  and 
Ad>N  =  10°.  The  results  showed  the  time-averaged  transmittance  to  be  equal  to  the  Fresnel  transmittance 
within  1  percent,,  f  1 10)  dt  =  t(0),  for  source  zenith  angles  0°to  65°.  The  values  for  zenith  angles  greater 
than  65°  are  given  in  Table  7.  The  slight  decrease  and  then  increase  with  wind  speed  at  zenith  angle  70° 
is  consistent  with  the  results  for  the  smaller  zenith  angles.  For  zenith  angles  0°to  60°,  the  time-averaged 
transmittance  consistently  decreased  with  increasing  wind  speed  Oto  38  knots,  even  though  the  decrease 
was  less  than  1  percent  between  0  and  38  knots.  The  values  in  Table  7  may  be  used  for  all  wavelengths 
in  the  visible  spectrum  fcr  both  salt  water  and  fresh  water  if  the  desired  accuracy  is  no  greater  than 
1  percent. 


Table  7.  Time-averaged  transmittance  for  n’  =  1.341  for  source  at  zenith  angle  theta. 


Zenith 

A°9le  Wind  Speed  (Knots) 

6 


Degrees 

0 

2 

4 

8 

14 

20 

26 

32 

38 

70 

0.864 

0.864 

0.863 

0.863 

0.865 

0.868 

0.872 

0.875 

0.879 

75 

0.785 

0.787 

0.790 

0.797 

0.809 

0.820 

0.830 

0.838 

0.845 

80 

0.649 

0.666 

0.680 

0.708 

0.739 

0.763 

0.781 

0.795 

0.807 

85 

0.414 

0.499 

0.542 

0.603 

0.660 

0.698 

0.726 

0.747 

0.765 

90 

0.0000 

0.359 

0.419 

0.506 

0.585 

0.636 

0.672 

0.699 

0.722 

3.7  Directional  Radiance  (Luminance)  of  the  Time-Averaged  Wind-rufflod  Sea 

The  directional  radiance  (luminance)  of  the  wind-blown  sea  due  to  time-  or  space-  averaged  small  ca¬ 
pillary  waves  is  the  sum  of  the  terms  for  sun  glitter,,  sky  reflectance,  and  upwelling  light.  In  radiometric 
form. 


Wp.  V  =  Ng(0p,  <f>p)  +  Nslty(0p)  +  Nu(0p,  c£p). 


(3.7.1) 


The  photometric  form  is  the  same. 


4.  RADIANCE  (LUMINANCE)  PATTERNS  OF  LARGE  WAVES 

Photographs  of  the  sea  surface  usually  contain  radiance  patterns  which  indicate  that  the  radiances  are 
not  time-  or  space-averaged.  The  photographic  work  by  Cox  and  Munk  (Ref.  15)  to  obtain  the  slope  distri¬ 
bution  utilized  a  defocused  camera  in  order  to  obtain  photographs  of  averaged  radiances.  The  resolution 


34 


of  many  photographs  is  such  that  the  radiance  pattern  indicates  the  larger  wave  train  pattern.  The  small 
capillary  waves  however  are  normally  not  resolved  and  can  be  assumed  to  be  time-  or  space-averaged.  The 
calculation  methods  described  in  the  previous  sections  can  be  easily  modified  to  give  a  reasonable  ap¬ 
proximation  of  the  intermediate  resolution  obtained  photographically.  This  is  done  by  tipping  the  ergodic 
cap  in  the  direction  of  the  normal  from  the  surface  of  the  large  wave.  The  surface  normal  used  should  rep¬ 
resent  the  average  slope  for  the  resolution  of  the  photograph,  such  as  1  ft.,  2  ft.,  4  ft.,  etc. 

The  equations  in  this  fourth  section  will  be  stated  in  radiometric  terms.  The  photometric  equations  are 
directly  equivalent. 


4.1  Tipped  Ergodic  Cap 

For  a  calculation  of  intermediate  resolution  of  the  time-  and/or  space-averaged  slope  distribution  of 
the  sea  surface,  it  is  desirable  to  have  the  capillary  or  small  facet  wave  distribution  separated  from  the 
larger,  gravity  wave  component.  The  capillary  distribution  probably  is  partially  dependent  upon  position 
on  the  wave  train,  at  or  near  the  peak  of  the  gravity  wave  the  distribution  would  probably  have  a  larger 
standard  deviation  than  down  in  the  trough.  No  such  statistics  being  available,  however,  it  is  probably 
not  unreasonable  to  assume  that  the  wave  slope  distribution  as  measured  by  Cox  and  Munk  is  primarily 
the  capillary  distribution  and  for  a  first  approximation,  can  be  assumed  to  be  the  same  for  all  positions  on 
the  larger  waves,  the  primary  effect  being  that  of  tipping  the  distribution  according  to  the  slope  of  the 
large  wave. 

Various  assumptions  may  be  made  about  the  point-to-point  slope  distribution  of  the  larger  waves.  The 
simplest  and  least  accurate  is  that  of  assuming  trochoidal  shaped  waves  in  the  direction  of  the  local  wind 
and  of  superimposing  trochoidal  shaped  swell  from  other  directions,  to  represent  swell  from  distant  storms. 
This  method  is  fully  described  in  Reference  3.  A  more  sophisticated  method  would  be  to  represent  the  sea 
surface  by  Lagrangian  equations.16 

In  tipping  the  ergodic  cap'in  the  direction  of  the  normal  from  the  larger  wave  surface,  the  vector  repre- 
senting  the  normal  from  the  wave  will  be  designated  as  T  with  the  zenith  angle  being  &T  and  the  azimuth 
from  the  sun  being  <f>T.  See  Fig.  22. 


i 


Kig.  22.  The  tilted  ergodic  cup. 
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The  vector  T  will  define  a  rotated  system  designated  by  double  primes.  The  normals  from  the  elements  of 
the  ergodic  cap  N  will  be  designated  by  angles  0“  and  0^  which  are  angles  in  the  double  prime  system 
[equivalent  to  angles  in  the  unrotated  system  when  0T,  <5 i>T  both  =  0,  when  the  cap  is  untipped].  The  hori¬ 
zontal  projected  area  of  the  element  representing  the  probability  of  occurrence  of  slope  is  in  the  double 
prime  system  and  is  designated  by  AH(0" ).  The  true  area  is 


At{0n  ) 


Ah(3n> 

cos0^ 


(4.1.1) 


4.2  Sun  Glitter 

The  sunlight  reflected  from  a  wave  tilted  at  0T,  <£T is  found  by  equations  similar  to  those  in  Section 
3.4.  Eq.  (3.4.1)  for  Ng(0p,  (f>p)  can  be  used  directly  but  the  equation  for  the  projected  portion  of  the  ergodic 
cap  reflecting  the  sun,  Ap(  At)  changes  from  Eq.  3.4.12  to 


tan’^ 

e 2  exp - 

2a2 

A  (At)  = -  (4.2.1) 

8a2  COS‘,0gN 


The  modified  equation  for  the  total  projected  area  of  the  cap,  Ap(  will  be  given  in  the  next  section. 

The  Fresnel  reflectance  angle  0SR  is  found  in  the  same  manner  as  for  an  untipped  cap  through  thc  use 
of  Eqs.  (3.4.14)  and  (3.4.16).  In  order  to  obtain  the  normal  from  the  cap  surface  reflecting  the  center  of 
the  sun  in  the  double  prime  system,  0g'N,  it  is  first  necessary  to  obtain  the  normal  position  in  the  unro¬ 
tated  system,,  #SN  and  cpSN.  0SN  can  again  be  found  through  the  use  of  Eqs.  (3.4.14)  through  (3.4.17).  In 
order  to  find  <£gN,  however,  it  will  be  necessary  to  also  compute  sin<^  by  appropriately  substituting  into 
Eq.  3.2.10.  As  before  the  prime  system  is  defined  by  the  vector  S.  Since  ^  =  c£g  =  0,  some  of  the  terms 
drop  out. 


sin^  =  [sin0osin<£j  -  sin0g 


(4.2.2) 


Then  <£SN  can  be  found  by  rotating  back  tc  the  principal  coordinate  system  and  appropriately  substituting 
into  Eqs.  3.2.12  and  3.2.13. 
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COS  (f> 


SN 


r  ■ 

=  sin  — 

L  2 


C0S<£o  COS0t  +  cos 


*“  1 

5 — sin0,  - 

2  J 


sin0, 


SN 


(4.2.3) 


r  e°  i 

sin$sN  =  [  sinysin^  j  -p 


sinft 


SN 


(4.2.4) 


Now  rotating  to  the  double  prime  system  defined  by  the  vector  T  at  0T,  <£T,  and  substituting  into  Eq.  3.2.8, 
f?s'N  is  obtained.  See  Fig.  23. 


cos0gN  =  sin(?SN  cos<£SM  sin0Tcos0T  +  sin#SN  sin<£SN  sin&rsin<£T  +  cos0SN  cos0T  (4.2.5) 


4.3  Sky  Reflectance 

The  sky  reflectance  for  the  tipped  wave  surface  is  a  direct  extension  of  the  equations  given  in  Sect. 
3.5  for  the  untipped  ergodic  cap  sky  reflectance.  Eq.  3.5.1  becomes 


n  r  m  "I  AH(0Ni) 

1  1  N(0b,  c6b), j  Q  UW  r(W  cos(W, 

■  ---  1  L  j  ~  1  J  cosf^'  t 


(4.3.1) 


The  expression  for  the  horizontal  area  of  the  surface  at  (%  =0  is  Eq.  3.5.2.  The  general  expression  for 
the  horizontal  area  of  each  element  merely  substitutes  6£  for  6^  in  Eq.  3.5.3.  Similarly  the  total  projected 
area  of  the  cap  changes  from  Eq.  3.5.5  to 


A 


P 


(4.3.2) 


The  Fresnel  reflectance  angle,  0BR,  and  the  sky  radiance  angles,  0B,  are  found  by  first  rotating  to 
the  double  prime  system  defined  by  the  tilt  position  0T.  <^..  The  observers  position,  0o,  <£o ,  is  found  in  the 
rotated  system  by  substituting  into  Eq.  3.2.8  through  3.2.10  as  follows. 


cos^'  =  sin0o  cos <j>o  sin0Tcos<£T  +  sin0o  sin<£o  sin0Tsin9ST  +  cos0o  cos0T  (4.3.3) 


cos<^'  =  [  sin0o  cos<£,  cos0Tcos0T  +  sin0o  sin0o  cos0Tsin^T  -cos0o  sin0T]  4-  sin6£'  (4.3.4) 


sin<^'  =  ( -  sin0o  cos<£o  sin<f>T  +  sin0o  sin0o  cos^>T]  *  sin6£' 


(4.3.5) 


See  Fig.  24  for  the  relationships  between  the  vector  0  and  the  unrotated  and  rotated  axes. 


Fig.  24.  Observer  O  vector  ,m  principal  and  double  prime  system. 
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Now  both  the  vector  0,  and  the  vectors,  N,  of  the  elements  of  the  ergodic  cap  are  in  the  double 
prime  system. 

The  second  step  is  to  let  the  vector  0  (at  0o,  <f>o  or  0o  ,  <f>0)  define  the  single  prime  system,  and  to  find 
the  vector  N  position  in  the  single  prime  system  through  substituting  into  Eqs.  3.2.8  through  3.2.10.  See 
Fig.  25. 


~  k 


■  l 

J 


Fig.  25.  Normal  N  from  element  of  ergodic  cop  in  double  prime  and  prime  coordinate  systems. 


cosf?^  =  sindj^  cos<f>x  sin6£  cos<^'  +  sinf?^  sin<£^  sin0°  sin<^'  +  cosO'^  cos6£'  (4.3.6) 

cos<£^  =  [  sin0"  cos cosd"  cos<^'  +  sin0^  sin^  cos^'  sin<£"  -  cosd^  sin£>£'  ]  *  sinf^  (4.3.7) 

sin<£^  =  [- sinO^cos^  sin<^'  +  sin^  sin^  cos^']  h-  sinO^  (4.3.8) 

As  before,  0^  is  the  Fresnel  reflectance  angle,,  0BR.  When  0'u  ^  90°,,  the  reflectance  and  projected  area 
are  zero  and  there  is  no  contribution  to  the  sum  from  that  element  of  the  ergodic  cap.  The  radiance  direc¬ 
tion  angles  in  the  prime  system  are  again  0B ■=  20^  and  <3^- 

The  third  step  is  to  go  from  the  single  prime  system  to  the  unrotated  system  through  the  use  of  Eqs. 
3.5.11  through  3.5.13  as  was  shown  in  Fig.  17,  section  3.5.  Again  when  0B^9O°,  the  radiance  contribution 
is  considered  negligible. 
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4.4  Upwelling  Radiancb 


It  will  be  assumed  that  since  the  untilted  ergodic  cap  represents  the  time-  or  space-averaged  case, 
that  it  is  also  appropriate  for  all  parts  of  the  upwelling  equation  relating  to  the  light  penetrating  the  sea 
surface  from  above,  and  reflecting  and  interreflecting  beneath  the  surface.  The  only  term  in  which  the  cap 
is  assumed  tipped  is  the  time-averaged  transmittance  of  the  radiance  through  the  water-air  surface  toward 
the  observer.  Thus  Eq.  3.6.1  becomes 


Nu«?p,  <f>p)  = 


£  ^  N(0b,  0b)Q °(0b) y“t(0B)cos0BdtdftB 
+  Ns(0,,O)Q°(0t) /t(0)cos6,df  dQ 

R«M')dt  /^y 

1  r,  dt)  \  n'  / 


1 

EJ 


rr  ( 1  -  R 


•y 


(4.4.1) 


The  angle  toward  the  observer  in  the  double  prime  system  0",  was  found  in  Eq.  (4.3.4)  in  the  pre¬ 
ceding  section. 


4.5  Directional  Radiance  of  Tilted  Sea  Surface 

The  directional  radiance  of  the  sea  surface  at  resolutions  which  do  not  average  over  large  areas  (or 
are  not  long  time  exposures),  can  be  approximated  by  computing  the  radiance  of  a  tipped  ergodic  cap.  Eq. 
3.7.1  is  directly  applicable  for  summing  the  components  of  the  radiance.  The  radiances  of  the  parts  of  the 
wave  pattern  would  emerge  from  computations  for  ergodic  caps  tipped  at  the  appropriate  angles. 
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APPENDIX  A 


DERIVATION  OF  EQUATIONS  FOR  THE  TRANSFORMATION  OF  ROTATION  OF 
AXES  EXPRESSED  IN  SPHERICAL  COORDINATES 

By  J.  W.  Wasserboehr  and  J.  I.  Gordon 


The  basic  equations  used  for  the  development  of  the  angular  relationships  were  those  for  the  transfor¬ 
mation  of  rotation*  of  axis  in  a  three  dimensional  rectangular  coordinate  vector  system. 

The  equations  are: 


x  =  x'(i  *i')  +  y‘(i  ‘j')  +  z'(i  -k') 

y  =  x'(j  • +  y* (J  *j’)  +  z'(j  ' k')  Group  (A.1) 

z  -  x'(k'i')  +  y'(k-j')  +  z'(k-k') 

x’  =  x(i  •  i')  +  y(  j  •  i')  +  z(k* i'J 

y'  =  x (i  -  j')  +  y(j  *  j')  +  z(k*  j')  Group  (A.2) 

z'  =  x (i  *k‘ )  +  y(j  *k')  +  z(k  -k') 


If  the  rectangular  coordinate  system  is  changed  to  the  spherical  coordinate  system  of  Figure  A.1  it  can  be 
seen  that. 


x  =  psin0cos?!> 
y  =  psin0sin<£ 
z  =  pcc,s0 


x'  =  psin0'  C0S9&' 
y '  =  psin0'  sin$’ 
z'  =  p  COS0' 


Group  (A.3) 


and  if  these  are  substituted  into  groups  (A.1)  and  (A.2),  and  the  common  factor  p  is  divided  out,  the  equa¬ 
tions  of  transformation  of  rotation  become. 


sin0cos<£  =  sin0'  cos<£'  (i  *i‘)  +  sin0'  sin<£‘  (i  ’j')  +  cos0'(i  *  k' ) 

sir«0  sin0  =  sin0'  cos.£'(i  ’i')  +  sin0'  sin<£‘  (j  -  j')  -1  cos0' (j  *  k' )  Group  (A.4) 

cos0  =  sin0'  cos (k  -  i')  +  sin0'  sin<£'  (k  •  j1)  +  cos0'  (k  •  k') 


*  J.  S.  Sokolni'ttoff,  and  E.  S.  Sokolnikoff,  Higher  Mathematics  for  Engineers  and  Physicists, 
McGraw-Hill,  N.  Y.,  2nd  Ed.,  1941,  p.  403. 


A? 


wppp, )  i  j  ,  n  m.u  1 1 1  u  - )  i  wjiiju.p '  i. » i  ii .! .  .p.w 


Fig.  A.  1.  Rectangular  Coordinate  System  to  Spherical  Coordinate  System. 
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sin0'cos<£'  =  sin0cos<£(i  •  i ‘ )  +  sin0cin<£(j  •  i')  +  cos0(k  •  i') 

sin0'  sin<£'  =  sin0cos<£(i  •  j')  +  sin0sin<£(j  •  j 1 )  +  cos0(k  •  j')  Group  (A.5) 

cos0'  =  sin0cos<£(i  ’k')  +  sin0sin<£(j  -k')  +  cos0(k  -k') 


In  spherical  coordinates  the  i,  j,  and  k  unit  vectors  are: 


i;  p  =  1.  0  =  90°,  </,  =  0° 
i;p  =  1,0  =  90°,  <f>  =90° 
k ;  p  =  1,0=  0 \<f>  =  0° 


From  these  it  is  obvious  that. 


i  =  (1)i  +  (0)j  +  (0,k  =  i 

j  =  (0)  i  +  (1)j  +  (0)k  —  j  Group  (A.6) 

k  =  (0) i  +  (0) j  +  ( 1 ) k  +  k 


If  0,  and  (£,  are  the  angles  for  the  transformation  of  rotation  to  the  new  axes  x',  y  _..d  z’  (Figure  A. 2) 
then  the  spherical  coordinates  of  i',  j',  and  k',  with  respect  to  the  original  x,  y,  and  z  systems  are: 
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i';  p  =  1.  6  =  0,  +  90°,  0  =  0, 

j';  P  =  1,  0  =  90°,  <£  =  0,  +  90° 

k' ;  p  =  1 ,  0  ~  ,  $  =  0i 


If  these  latter  values  are  substituted  into  the  general  equation  of  a  vector. 


A  =  xi  +  yj  +  zk  =  (psin0cos$)i  +  (p  sin0  sin0)  j  +  (pcos0)k 


then  equations  for  the  vectors  i ' .  j‘,  and  k1  in  relation  to  the  original 


system  are  obtained  as  follows: 

i'  =  [sin(0,  +  90)cos<£t]i  +  [sin  (0,  +  90)  sin0,]  j  +  [cos(0,  +  90)]k 

j'  =  [sin  90  +  cos(<£,  +  90)] i  +  [sin  90  sin (0,  +  90)]  j  +  [cos  90] k 

k'  =  [sin0, , cos0,]  i  +  [sin@,  sin0,]  j  +  [cosf?,]k 

and  sin.  sin  90°  =  1,  cos  90°  =  0,  sin(0,  +  90)  =  cos0,;  and  cos(0,  +  90)  =  -sintf, 

i'  =  (costf,  cos<£,)i  +  (cosf?,  srn0, )  j  -  (sin0,)k 

j'  =  -  (stn0, )  i  +  (cos0,)j  +  (0)  k 

k'  =  (sin0,  cos0, )  i  +  (sinft,  sin0,)j  +  (cosfl,)k 


Group  (A.7) 


Group (A. 8) 


Group  (A. 9) 
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From  vector  analysis  the  dot  product  of  any  two  vectors  is  A  •  B  =  xa  xb  +  yayb  +  za  zb,  therefore, 

i  *  i '  =  (1)  (cos6lcos<f>l)  +0  +  0  =  cost?,  cos0, 
i  *  j'  =  (1)  (-sin^j)  +0+0=-sin^1 

i  *k'  =  (1)  (sin^cos^,)  +  0  +  0  =  sin^cos^, 

j  *  j'  =  0  +  (1)  (cosfl,  sin0,)  +  0  =  costfj  sin0, 

j  •  j  *  =  0  +  (1)  (cos^,)  +0  =  003$!  Group  (A.10) 

j  *k‘  =  0  +  (1)  (sin0,  +  sin0,)  +  0  =  sin0,  sin0, 

k  *  i'  =  0  +  0  4  (1)  (-sin^)  =  -sin#! 
k  *j‘  =  0  +  0  +  (1)(0)  =  o 

k  *k'  =  0  +  0  +  (1)  (cos0j)  =  cos0, 

Substituting  these  dot  product  values  into  the  group  equations  (A.4)  and  (A.5),  the  following  are  obtained. 

sin0cos0  =  sinfl'cos^'costfjcos^,  -  sin0'  sin0'  sin0,  +  cos0'  sin0j  cos^ 
sin0sin0  =  sin0‘ cos0' cos0,  sin0,  +  sin0' sin0’ COS0,  +  cos0' sin0j  sin0,  Group  (A.11) 

cos 0  -  -sin0'cos0' sin0j  +  cos0'cos0, 

and, 

sin0'cos0’  =  S!n0cos0cos0!cos0j  +  sin0sin0cos0,  sin0,  -  cos0sin0, 
sin0'sin0'  =  -sintfcos^sin^  +  sin0sin0cos0j  Group(A.12) 

cos  O'  -  sin0cos0sin0!  COS0J  +  sin0sin0sin0,  sin0,  +  cos0cos0( 
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APPENDIX  B 


OCCLUSION 

by  Keith  B.  MacAdam 


The  theory  of  occlusion  is  given  in  Reference  2,  Chapter  12.  The  problem  will  be  stated,  definitions 
and  terminology  given,  and  some  results  derived. 

Definition 

Q°(x,  £)  is  the  average  fractional  time  that  an  incident  path  £  corresponds  to  direct  light  from  the  sky 
or  sun  at  point  xon  the  mean  sea  surface  &  As  the  waves  undulate,  the  intersection  x  of  the  path  of 
sight  £with  S  changes  its  horizontal  position;  Fig.  B.l. 


Fig.  B.l 


On  the  assumption  that  the  sea  surface  elevations  C  at  neighboring  points  x  are  uncorrelated,  it  can 
be  imagined  that  what  is  seen  from  f  is  the  air-water  interface  point  directly  above  or  below  x,  with  no 
change  in  results.  Fig.  B.2. 
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Fig.  B.2 


The  problem  can  be  visualized  as  follows:  How  often  could  one  see  a  floating  bob  at  horizontal  point 
x?  (Indeed,  this  point  of  view  yields  an  immediate  technique  for  empirical  verification  of  the  theory). 

There  is  a  certain  probability  pt(£)A£  that  x  will  fall  into  the  vertical  interval  (£,£  +  A£).  Wher¬ 
ever  it  falls,  there  is  a  certain  probability  that  £has  n'o  additional  intersections  with  S  between  x  and  the 
sky. 


If  p2(£,  A£)  is  the  probability  that  in  an  increment  A£  of  vertical  ascent  along  £  there  is  no  inter¬ 
section,  then 


IlPgK.AO 


is  the  probability  that  f  is  completely  free  of  intersections.  Then  the  final  answer  is 


Q°  =  2Pl(£)A<  [llPzUr.AO] 


(B.1) 
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In  the  following,  it  is  assumed  that  the  Neumann  spectrum,  developed  in  Ref.  2  satisfactorily  repre¬ 
sents  the  ocean  surface,  although  it  is  essential  only  that  the  mathematical  form  of  the  Neumann  spectrum 
agree  with  experiment. 

From  the  Neumann  spectrum. 


P,(£)A£ 


(B.2) 


where  m  is  the  mean  square  wave  height. 


From  the  Neumann  spectrum  it  is  also  possible  to  obtain  the  average  number  of  times  per  unit  length 
that  a  horizontal  path  at  height  £  from  the  mean  sea  surface  crosses  the  air-water  interface.  That  is, 


n(£) 


i  r  m2o  "l 

Lm.0  J 


exp 


_£^ 

2tn 


(B.3) 


where  m20  =  a2,  the  mean  square  wave  slope.  Assuming  that  ideas  from  the  Poisson  distribution  describ¬ 
ing  random  events  may  be  usud,  the  probability  that  a  horizontal  path  of  length  Ax  is  free  of  intersections 
is  exp  (-n(£)  Ax). 

Then  if  the  slanted  path  £  is  approximated  by  a  series  of  short  horizontal  paths  at  increasing  height, 
the  result  is  the  approximation 


riP2<£.A£)  =  n  exp  A£/tan</f3  (B-4> 

j  =  o 


since  A£  =  Ax  tan^.  4*  is  the  elevation  angle.  This  is  a  procedure  similar  to  that  used  in  Ref.  2.  Then 


IIp2<£.A£) 


n  exp  [-n(£o  +  jAxtan</<)  A£/tani/»] 

j=0 


(B.5) 


=  exp  \  1  ~n(£0  +  jAxtaniA)  A£/tar.i/> 
f  j  =  o 


(B.6) 


Now  setting  £  =  £o  +  j  Ax  tarn//,  and  letting  Ax-»0 


n=  exp 


(B.7) 
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Upon  evaluation 


49 


r 


<V 


\ 


(9.16) 


Then 


(B.17) 


the  effect  of  the  extra  factor?  For  suns  higher  than  «A  -  mte9ra'  With  va'Ue  h  How  lar9e  is 

integral  is  less  than  12%  greater  than  1  and  Drobablv  h  »  'J”  ^  W°rSt  C3Se'  V  =  38  knots'  the  actual 

than  1.08.  For  suns  higher  than  0  =  20°  the  inteoral  6  7  1  **  V  =  26  th®  3CtUal  integral  is  ,ess 

equal  to  38.  *  **  13  less  than  103  «H  wind-speeds  less  than  or 


Then  for  10°  <  ^  <  90°  and  0  <  v  <  38  the  formula  below  may  be  used 


Q°  =  exp  [-o/(/ Yn  tan 
It  will  be  good  to  within  12%.  In  all  cases  the 


(8.1 8) 


within  the  specified  lir.iit. 


actual  value  of  Q°  will  be  greater  than  that  given  above. 


C  S  ,led  ,0  ,he  s“rf“e  rather  than  the  actaal  surface  8.18  It 


tamed  exactly. 


is  ob- 
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